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Abstract 

We analyze the decays (jj KK utilizing a formulation of transition rates which explicitly exhibits 
corrections to Fermi’s Golden Rule. These corrections arise in systems in which the phase space and/or 
matrix element varies rapidly with energy, as happens in (/ —> KK, which is just above threshold. We 
show that the theoretical corrections resolve a puzzling 5a discrepancy between theory and experiment for 
the branching ratio R = r((/ ^ A+A“)/r((/ ^ K^K^). 


One of the most well known results from elementary quantum mechanics is the formula relating the rate T for the 
transition \i \n > to the matrix element Vni{E) =< n\V\i > induced by a time-independent perturbation V, 

Ott 

r(2^n) = —lV„.(E)l^p(E), ( 1 ) 

where E = En — Ei, and p{E) = dN/dE is the energy density of the final states n [1]. This formula is so useful and so 
widely applied that Fermi named it “Golden Rule No. 2” (FGR2) [2], As almost all derivations of Eq.(l) make clear, 
FGR2 is an approximation which is valid in the limit when both \ Vni{E)\‘^ and p(i?) are “slowly varying” functions of A, 
although the precise meaning of “slowly varying” is not always made explicit. However, it is evident that \ Vni{E)\‘^p{E) 
will not be slowly varying in some circumstances, for example, when the initial state \i > is just above the threshold 
for decay into the final state |n >. This is the case for the decays </> ^ KK (either (p K~^K~ or p ^ K^K^), where 
= 1019.417(14)Mey, while TOk+ = 493.677(16)MeR, and niKO = 497.672(31)MeR [3]. Quark model diagrams 
for these decays are shown in Fig. 1. In principle the corrections to FGR2 in such decays could be large enough to 
lead to detectable effects, and in what follows we show that this is in fact the case. More interestingly, we demonstrate 
explicitly that the correction to FGR2 arising from the rapid variation of \Vni{E)\^p{E) with E resolves a puzzling 
discrepancy [ 5 ] between theory and experiment for the ratio R = r{p K~^K~)/r{p 77 ° AT^). 

An explicit expression for the corrections to FGR2 can be conveniently derived by endowing the initial decaying 
state at the outset with a lifetime r = l/F, and then solving self-consistently for F. Consider an initial state \i,to > 
which evolves into the state \i,t > at a later time t under the influence of the time evolution operator U{t,to). The 
state \i,t> can be expanded in terms of a complete set of eigenfunctions |n > of the unperturbed Hamiltonian, 


\i,t> = |n >< n\U{t,to)\i,to >= y^ c„{t)\n >, 

n n 

Cn(t) = ^ [ 

" Jto 

where huni = En — Ei, and Ei{En) is the unperturbed initial (final) energy. Without loss of generality we can set 
to = 0, the instant at which the p is produced. The quantity of interest is c„(oo) which is given by 


C„(oo) = -^ / dt'e^"‘^''' Vni = l-Vr. 


n ”Qw™-r/2’ 


( 3 ) 


where Vni is independent of time. We next impose the unitarity constraint on Cn(oo), namely |c„(oo)p = 1. After 
the sum is converted into an integral in the usual manner, ^ f dEp{E), the unitarity constraint assumes the 
form 


1 = 



dEp{E) 


\V{E)\^ 

A2 -h (?ir/2)2 


( 4 ) 


where we have set fiLOni E, and Vni E(i?). The denominator in Eq.(4) is rapidly varying in the vicinity 
of A = i?o = 0, which corresponds to an energy-conserving transition. Thus if we invoke the assumption that 
p{E)\V(E)\^ is slowly varying with respect to the denominator [E'^ + (?ir/2)^], then the unitarity constraint in Eq.(4) 
yields 

l^piEo)\V{Eo)\^ dE^^^j^^ = p{Eo)\V{Eo)\^ ■ i2n/hr). (5) 


Solving Eq.(5) for F we are led immediately to the standard Fermi Golden Rule in Eq.(l). Moreover, the Golden 
Rule integral technique (GRIT) embodied in Eq.(4) gives a specific formula for the corrections to FGR2 in Eq.(I) for 
processes in which p{E) and/or |E(if)p varies significantly with energy. One can further elucidate the approximation 
being made in going from Eq.(4) to Eq.(5) by invoking the identity 


lim 


0 f;2 




= -m, 

a 


( 6 ) 


where a = hT/2. Gombining Eqs.(4) and (6) leads immediately to Eq.(5) in the limit F ^ 0. However, for F 0 
5{E) is replaced by the (broader) Lorentzian in Eq.(4) which introduces contributions from p{E)\V{E)\'^ in which 
A yf 0. As we now demonstrate, these additional contributions are quantitatively different for p K~^K~ and 


2 



(j) K^K^, and their inclusion serves to resolve a 5a discrepancy between the theoretical and experimental values 

[4,5] of i? = r((/) ^ K+K-)/ri(j) K°K°). 

The decay </> ^ K^K~ is induced by the Lagrangian density 

£(x) = ig+(j)^{x)[K+ {x)d^j,K~{x) - K~{x)d^K+ {x)], (7) 

where is the appropriate coupling constant, and K^{x) annihilates , etc. A similar expression characterizes 
(j) —> which is proportional to the coupling constant go, with go = g+ in the limit of exact SU(2) symmetry. 

Bramon, et al. [5] have considered the effects of radiative corrections, and we will return to their results below. Using 
FGR2 as given in Eq.(l) the decay rate r{(f) K^K~) obtained from Eq.(7) is given by (setting h = c=l hereafter), 

rw - k * k -) = H (I) 1|, (8) 


where jfcj = (1/2 )(to^ — is the magnitude of the 3-momentum in the </> rest frame. The factor of can 

be understood as follows: Since the kaons are spinless, whereas (() is a vector particle, angular momentum conservation 
demands that K and K be emitted in a relative P—wave, which is consistent with the derivative coupling in Eq.(7). 
Hence |< Vni >P contributes a factor of |fc|^, while the phase space contributes an additional factor dN/dE oc |fc|. 
Combining Eq.(8) with the corresponding expression for (j) we find using FGR2, 


_ T{(j>^K+K-) 



RfGR2 


(gl\ / I -44 

UoV Vl-4/i^ 


3/2 


(9) 


where /i+ = mj^+/m^ and /ig = nixo/m^. Inserting the previously quoted values of mx+,rnj^o, and m,p into Eq.(9), 
and assuming go = g+, we find Rth = 1.528, to be compared with the experimental value [4] 


Pexp = 1.456 ±0.033 . (10) 

Bramon, et al. [5] have evaluated various corrections to Rth in an effort to bring Rth and Rexp into agreement. 
Most significant among these are electromagnetic radiative corrections which affect </> ^ K~^K~ but not (p . 

These authors find that the radiative correction factor, 77 = 1.042, increases Rth to 1.59. Bramon et al. have also 
studied the effects of SU(2) symmetry breaking on the ratio g+jgo, which arise via quark mass differences. The 
(j) wavefunction is pure ss, and hence the K^K~ (K^K^) final state requires the creation of an additional uu(dd) 
pair (see Fig. 1). Since the uu pair is lighter than dd, one expects this effect to enhance </> ^ K^K~ relative to 
(j) , thus further widening the discrepancy between theory and experiment. A detailed analysis by Bramon 

et al. [5] finds g+jgo — 1-01, which agrees with the intuitive expectation that this correction also works in the wrong 
direction. With this correction included Rth becomes 1.62. Other effects considered by these authors, such as the 
inclusion of electromagnetic form factors in calculating radiative corrections, and final-state rescattering effects, are 
negligible. We are thus left with a puzzling 5cr discrepancy between Rth = 1-62 and Rexp = 1.456 (33). 

We proceed to demonstrate that this discrepancy can be resolved by incorporating the corrections to FGR2 that 
arise from the Golden Rule integral technique. Gombining Eqs.(4) and ( 8 ), and introducing the notation 2 ; = E/m^, 
a = VI2m^, (F = 4.458(32)MeV) we express the unitarity constraint for (/(-decays in the form, 


1 


1 

37r 



dz 1 

( 1 ± 2 ) z2+a2 





1 

37r 



dz 1 
{l + z) z'^ + a? 




-I 3/2 


( 11 ) 


The two terms exhibited in Eq.(ll) are, respectively, the contributions from the K~^K~ and states, and ... 

denotes contributions to the unitarity integral from other channels (such as pir) which can be ignored for present 
purposes. We emphasize that the functional form of the expressions in square-brackets in Eq.(ll) is determined by 
the kinematics of (/> ^ KK, specifically by the relation between |fc| and E given in Eq.(12) below. The density of final 
states is readily found to be proportional to (1 ± z)k, and each of the bosonic normalization coefficients {2Ek)~^^'^ 
contributes a factor (1 ± 2 :)“^, resulting in an overall factor (1 ± 2 ;)“^. Eq.(12) also fixes the lower limit of integrations 
in Eq.(ll) as we discuss below. Rth is given by the ratio of the two terms in (11), which in the narrow resonance 
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(a 0) limit gives the standard result in Eq.(9). To specify the integration limits we note from Eq.(4) that since 
E = hojni is the energy difference between the final state \n > and the initial state \i >, we can write for (j) K'^K~ 

in the (j) rest frame, 

E = 2y |fcp + - TO0. (12) 

The lower limit on E evidently corresponds to |A:| = 0, and gives Emin = ‘2rnx+ — rn^. Accordingly in Eq.(ll), 
Zmin = (2/r+ — 1) for (j) K'^K~, and Zmin = (2^0 “ 1) for (j) The upper limit on |A:| (and hence z) 

extends to infinity. This limit leads to divergent integrals in Eq. (11), so the unitarity constraint (here as elsewhere) 
requires modification of the high-energy behavior of the <j)KK amplitudes. This can be achieved by incorporating a 
phenomenological form factor. 


F{W) 


M2 -h |fc|2’ 


(13) 


multiplying the (j)KK amplitude. This form factor introduces an asymptotic 1/z^ dependence (after the <j>KK 
amplitude is squared); so convergence is assured. The energy scale, M, is related to the confinement size of the 
hadrons involved (compared to 1/fc), and is typically of order ~ 1 GeV [6]. We have calculated Rgrit numerically as 
a function of M, and combined those results with the radiative correction factor rj = 1.042 and the SU(2) correction 
{ 9 +/ 90 ) = 1-02 to obtain Rth, 


Rth = ^vRgrit ■ (14) 

9o 

A plot of Rth as a function of M is shown in Eig. 2, along with the Icr experimental result from Eq.(lO) which is 
indicated by the dashed horizontal line. We see from this figure that Rth is relatively insensitive to the choice of M, 
and that for M ~ 0.8 GeV Rth falls within the Icr experimental bounds. Eor the nominal value M = 1 GeV we find 
Rth = 1.48 compared to the experimental value Rexp = 1.456 (33). 

The reduction in R relative to its value derived from EGR2 can be understood by considering the integrand of the 
integral in Eq.(ll), shown in Fig. 3. The factor multiplying the Lorentzian denominator is asymmetric about 
z = 0. The contribution from z > 0 significantly exceeds the result obtained if this factor is replaced by its 2 = 0 
value. The proportionate increase is greater for the decay than for the K^K~ decay because </> ^ is 

closer to threshold, so R becomes smaller than the value in Eq.(9). 

Although the discrepancy between the theoretical values of Rgrit and Rfgr2 is ~ 9%, the corrections to the 
individual partial decay rates are larger. Tgrit/^fgr2 is shown as a function of M in Fig. 4 for both K^K~ and 
channels. It seems evident that corrections to FGR2, similar to those considered here (but not necessarily so 
dramatic) can be anticipated in other decays. 
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FIG. 2. The solid curve gives the theoretical branching ratio versus M, the energy scale in Eq.(13), where 
Rth = {g+ / goYrjR-GRiT ~ I.QQZRgrit ■ The experimental ratio in Eq.(lO) is given by the long-dashed line, and the Icr 
bounds fall within the dotted lines. The short-dashed line is the result predicted by FGR2. 
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